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Abstract. We review models of cosmological gamma-ray bursts (GRBs). The statistical 
and 7-ray transparency issues are summarized. Neutron-star and black-hole merger sce- 
narios are described and estimates of merger rates are summarized. We review the simple 
^\ fireball models for GRBs and the recent work on non-simple fireballs. Alternative cos- 
Q\ mological models, including models where GRBs are analogs of active galactic nuclei and 
r- H where they are produced by high-field, short period pulsars, are also mentioned. The value 
of neutrino astronomy to solve the GRB puzzle is briefly reviewed. 
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ON 1. Statistical Premonitions 

[- — Soon after the first published report (Klebesadel et al., 1973) of the discov- 
ery of gamma-ray bursts (GRBs), Usov and Chibisov (1975) discussed the 
statistical issues. The size distribution N(> F) for GRBs produced by galac- 
tic disk neutron stars (NSs) bends from N(> F) oc F -1 for GRB fmences 

IT) F(ergs cm" 2 ) < F b to N(> F) oc _F" 3/2 for F > F b , with strong accompa- 
nying anisotropy of the GRB sky distribution towards the galactic disk when 

^ F <C F b . For burst energies Q = 10 38 <538 ergs and scale heights of 300g? 2 i 
pc, the break F b in the fmence size distribution occurs at F b = 10 _5 <538^2i 2 - 
O Two cosmological models were considered by Usov and Chibisov, and the 

£j point was made that a bend in the size distribution due to cosmological 
redshifting would occur for sources located at redshift z = 1, or source 
distances d(cm) = c/Hq ~ 10 7 c?2i = c?28- For GRBs from supernovae with 
Q 38 = 10 10 , a flattened size distribution N(> F) oc k, a constant, at low 
fmences followed solely from expansion effects. But with the bend in N(> F) 
found at fmences F < F b ~ 10~ 9 ergs cm -2 , it would be be quite difficult 
to detect due to the background galactic 7-ray glow. On the other hand, 
if Q ~ 10 52 ergs, then F b ~ 10~ 5 ergs cm -2 . The proposed model for this 
incredible energy release was the collapse of magnetic supermassive stars 
with masses M ~ 1O 5 -1O 6 M0, which was then being considered (Ozernoi 
and Usov, 1973; Prilutskii and Usov, 1975) as a model for the engine of 
active galactic nuclei (AGNs). 

In 1983, van den Bergh published an early (the first?) Aitoff projection 
of the sky distribution of 46 GRBs. From their apparent isotropy, he con- 
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eluded that the bursting sources were either non-exotic galactic black holes 
(BHs) or NSs with limiting detector sampling distances dj ~ 0.1 — 1 kpc, or 
cosmological objects with z £ 0.1 in order to be in accord with the absence 
of GRB associations with M31, M33, SMC, LMC, Virgo, and other nearby 
galaxies. He therefore argued that dj £ 10 27 cm for cosmological sources. 

2. Fallout from the BATSE Discovery 

The initial BATSE result (Meegan et al., 1992) showed strong inhomogeneity 
to limiting peak count rates C p = 1 ph cm -2 s _1 and limiting peak fluxes 
(f) p = 10~ 7 ergs cm -2 s _1 in the 50-300 keV band. In the discovery paper, 
{V/Vm ax ) = 0.348 ± 0.024. As reported by M. Briggs at this conference, the 
analysis of 657 GRBs gives (V/V max ) = 0.330 ± 0.011. Briggs also reports 
that the dipole moment (cos 9) deviates less than +0.9<r from 0, and the 
quadrupole moment (sin 2 b) deviates by less than —0.3<t from the value 1/3. 
When the strong inhomogeneity of sources is coupled with burst directions 
in statistical accord with isotropy, galactic disk NS models are instantly 
demolished. 

High-velocity NS models (e.g., Li and Dermer, 1992; Lamb, 1995; T. 
Bulik, this conference) remain an option, but require fine tuning of 3 param- 
eters: (i) the sampling distance or intrinsic GRB luminosity; (ii) a delayed 
turn-on parameter; and (iii) strong suppression of the bursting rates of slow 
(v ^ 800 km s _1 ) NSs. Tuning of parameter (iii) is particularly improbable, 
so that the statistical evidence for isotropy favors the cosmological model 
by wide margins. 

3. Statistics Post-BATSE 

But do the simplest cosmological models fit the observed size distribution of 
GRBs? In terms of the peak flux size distribution in the energy range 100- 
500 keV (see Fenimore et al, 1993a), N(> (j> p ) oc </>-°- 8 for 10" 7 < (j> p (ergs 

cm -2 s _1 ) < 2xl0 -6 , and N(> (f) p ) oc (f) p 1 - 5 ( ±0 - 1 ) a t larger values of (f) p when 
combined with the PVO data. The earliest analyses of the statistics of GRBs 
after the announcement of the original BATSE results assumed no source 
evolution with z for either luminosity or comoving number density, mono- 
luminous GRBs, power-law or broken power-law spectra, and unbeamed 
sources that uniformly emit in all directions. Good fits to the BATSE data 
were in all cases possible, with model fits implying limiting redshifts zb 
of the faintest detectable BATSE GRBs, peak GRB luminosities L p and 
bursting rates p. Table 1 gives the model results for the early analyses. The 
main difference between these models is in the implied values of L p , which 
depends on the assumed bandwidth. The high value of Dermer (1992) is due 
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to the use of a broken power-law spectrum and a GRB bandwidth extending 
to 100 MeV, whereas the other analyses restrict L p to the BATSE bandpass. 

Wickramasinghe et al. (1993) considered different cosmological models 
with Q = 0.1 and Q = 1.0, but could not discriminate between the two cos- 
mologies from the BATSE GRB data. Tamblyn and Melia (1993) performed 
a K-correction for different GRB detectors, their associated bandwidths, 
limiting thresholds, and trigger criteria. A broken power law GRB spec- 
trum was used to approximate the curved GRB spectra, and reconcile the 
detection rates by PVO, SMM, KONUS, SIGNE, APEX, and BATSE. 

TABLE I 



Inferred Properties of Mono-luminous, Non-Evolving, Unbeamed GRBs. 



Reference 


ZB 


L p (ergs 


s" 1 ) 


Bursting rate p 


Piran (1992) 


1 


10 50 (1 s 


GRB) 


3 x 10" 7 yr" 1 galaxy" 1 


Mao and Paczyriski (1992) 


1.5 


2 x 10 51 




2 x 10" 6 yr" 1 (i, galaxy)" 1 


Dermer (1992) 


1.2 


4 x 10 51 




10" 7 yr" 1 Mpc" 3 


Fenimore et al. (1993a) 


0.8 


5 x 10 50 




2.4 x 10" 8 yr" 1 Mpc" 3 



Burst luminosity functions provide too much parameter freedom, and 
are therefore not well constrained. The only evident restriction is that the 
luminosity function be not sufficiently broad and gently varying so as to 
erase the rather abrupt change in slope in the observed size distribution 
from a —0.8 slope to a —1.5 slope. Beaming effects have been considered 
by, for example, Yi (1993) and Yi and Mao (1994). The flux density S (ergs 
cm" 2 s" 1 MeV" 1 ) observed from a plasma blob travelling with Lorentz 
factor r and speed /3c = (1 — r~ 2 ) _1 / 2 c with respect to the stationary frame 
in the Hubble flow goes as S oc L com V 3+a . Here it is assumed that the 
source isotropically radiates a spectrum with energy spectral index a and 
total photon luminosity L com in the comoving frame. The Doppler factor 
V = [r(l — /3/x)]" 1 , where 9 = arccos /x is the angle between the jet axis and 
the line-of-sight to the observer. Values of (V/V max ) ~ 0.33 are found when 
L com = 10 51 /r 3+a ergs s" 1 , which corresponds to fairly weak luminosities 
indeed if T ~ 10 2 — 10 3 . The total number of sources must be increased by 
a factor ~ 4T 2 to take into account those which are not directed into our 
viewing direction. 

4. Getting the 7 Rays Out: The Need for Speed 

Table 1 shows that at distances d = 10 28 c?28 cm > GRBs must have source 
luminosities 
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£(ergs s" 1 ) = d 2 (j>8{l « lO^C^-e^ , (1) 

where £ is a bandwidth correction factor, (j> = 10 _6 ^>_6 ergs cm -2 s _1 is the 
measured energy flux, and SQ is the solid angle into which the emission is 
beamed. From the Elliot- Shapiro (1974) relation, updated for Klein-Nishina 
corrections on the radiation force (Dermer and Gehrels, 1995), this immedi- 
ately rules out unbeamed Eddington-limited accretion models by ~ 9 orders 
of magnitude for a variability time scale St v = 10 ms. 

Gamma rays can only escape from the source region if its pair production 
optical depth r 77 through the reaction 77 — >e + e~ is <C 1. For MeV 7 rays, 
this is essentially equivalent to requiring that the photon compactness for 
unbeamed sources 

L a T _ 10 12 C4^- 6 < 1 . m 
i?m e c 3 ~ St v (s) ~ ' K} 

hence the formation of a pair fireball is inevitable for stationary sources at 
cosmological distances (Piran and Shemi, 1993). 

Bulk relativistic motion relieves 7-ray opaqueness in two ways (Krolik 
and Pier, 1991). For plasma blobs radiating isotropically in the comoving 
frame, the comoving frame photon energy e com (= /iz/ com /m e c 2 ) is reduced 
by a factor T from the observed energy e Q b s , and the inferred photon energy 
densities in the comoving frame are reduced by a factor r 3+a compared to 
inferences from stationary sources. To avoid 77 attenuation, it is sufficient 
to require that e com < 1 to avoid the pair production threshold, so that one 
requires that T > the maximum observed photon energy e b s ,max- In view of 
high energy photons observed with SMM and COMPTEL, and even > GeV 
photons with EGRET (e.g., Hurley et al. 1994), it seems more reasonable 
to consider power-law photon spectra extending without high-energy limits 
in the source. Baring (1993) shows that for selected GRBs, the reduction in 
the comoving frame energy density implies that T must be £ 10 2 to avoid 
77 absorption. 

One must also make sure that photons emitted at earlier times do not con- 
stitute an additional source of opacity. Fenimore et al. (1993b) and Woods 
and Loeb (1995) have considered an expanding shell geometry and find again 
that values of T ~ 10 2 or greater are again needed to get the 7 rays out. 
Everything points to expansion at relativistic speed. 

5. Simple Fireballs 

Cavallo and Rees had already treated in 1978 a system which naturally 
produced plasma expanding relativistically outwards. The most important 
processes in their cosmic fireball were pair production and annihilation and 
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Compton scattering. They identified an important quantity that related the 
total fireball energy Qf to the baryon mass M through the expression 

Q f =riMc 2 . (3) 

The term 77 is essentially the entropy or photon energy per baryon, or the 
baryon-loading parameter. Because the important parameters in the fireball 
were a baryon or proton optical depth r p = n p (jyiZ , and t 77 or the fireball 
compactness I oc L/Rqc, they identified the four fireball types shown in 
Table 2. The term Rq refers to the radius of the region into which Qf is 
injected, and n p is the mean proton number density. 



TABLE II 

Simple Fireball Types (Cavallo and Rees, 1978). 



r 77 > 1, t p > 1 
Baryon-Dominated Fireball 

r 77 < 1, t p > 1 
Thermalized Radiation 



r 77 > 1, T p < 1 
Pure Fireball 

r 77 < 1, t p < 1 
Line-Production Region 



For the baryon-dominated fireball, the terminal speeds {3 t c or terminal 
Lorentz factors T t of an expanding fireball were found to be 

fa * r? 1 / 2 , 77 < 1 (4a) 

r t ^ 77 , 77 >1. (46) 

Shemi and Piran (1990) introduced the radius -Rthin where the fireball becomes 
optically thin, so that Eq. (4b) is generalized to 



_|_ ^ ( -Rthin/-Roj if V > -Rthin/-Ro 

«i 77 [l + V, if V < Rthin/Ro 



(5) 



Although Cavallo and Rees had made the connection of fireballs to GRBs, 
it was Paczyhski (1986) who pointed out two remarkable coincidences between 
fireballs and GRBs. The first was that if some small fraction ~ 0.1% — 1% of 
the 10 53 — 10 54 ergs available as the rest mass energy of a NS were liberated 
as 7 rays from cosmological distances, then a fluence F ~ 10 _6 <55i<i^ 8 2 ergs 
cm -2 would be measured near Earth. NS phase transitions and coalescing 
NSs were proposed as possible burst origins. The second coincidence was 
that if this energy were injected into a source size R ~ 10 km on a time 
scale Sti of a second or less, then the blackbody temperature 
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T„(MeV),3 ff;/;^:; 2 , (6) 

nicely in the 7-ray range. 

Paczyhski solved the steady-state spherically symmetric flow equations 
for continuous injection of energy with no baryon loading. For adiabatic 
expansion, T(R) = R/Ro, and the comoving frame temperature T com (R) = 
TqRqI R. An observer would see Doppler-shifted emission with effective tem- 
perature T Q bs — T(R)T com (R) ~ To, so that the emission from the expanding 
fireball would also peak at 7-ray energies. Paczyhski recognized that a ther- 
mal spectrum does not give good fits to GRB spectra, and suggested that 
a cool envelope would scatter 7 rays to lower energies to produce a photon 
flux c/)(e) 25 e _1 , as observed. 

Goodman (1986) considered a more realistic time-dependent bursting 
source. The calculated spectrum is slightly broader than a Planckian. He also 
pointed out that the minimum variability time scale would be ~ AR/c ~ 
Ro/c and that the GRB duration would be ~ iZthin/TjC, which could in 
principle be as short as Rq/c. 



6. Merger Rates and Scenarios 

The discovery of the binary radio pulsar system PSR 1913+16 (Hulse and 
Taylor, 1975) with a merger time scale t meTge ~ 3 X 10 8 yr <C the Hubble 
time tn ~ 10 10 yrs, prompted Clark et al. (1979) to make a first estimate of 
NS-NS merging rates. Because there was one such system for the 300 then- 
known isolated radio pulsars, and because the NS birthrate ./Vns ~ 0.02 
yr -1 , this implies a birthrate of binary NSs equal to ./Vns-ns ~ 10~ 4 yr -1 . 
This is also the merging rate in steady-state if t meTge <C £h- 

Narayan et al. (1991) performed a better treatment of selection biases 
in binary and isolated pulsar searches, and also considered the additional 
information provided by 3 new binary pulsars. They estimated TVns-ns ~ 
10 _5 /io(kpc) yr -1 , where the unknown binary pulsar scale height ho is esti- 
mated at a few kpc. From the observed massive X-ray BH binary systems 
such as Cyg X-l and LMC X-3, they also estimated the birthrate for NS- 
BH binaries to be ./Vns-bh ~ 10~ 4 ' 5 yr -1 . For 10% of the binary NS sources 
having t meTge < tn, this implies -/V™g^ s ~ 10~ 6 yr _1 , which is consistent 
with the implied GRB burst rates in Table 1 provided that the sources are 
not strongly beamed. 

Tutukov and Yungelson (1993) estimate higher merger rates than Narayan 
et al. (1991). For our galaxy, they find that iV™ g e I^ s f» 3 X 10~ 4 yr -1 and 
^NSS-BH ~ 10~ 5 yr -1 . This shows that order-of- magnitude uncertainty still 
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remains. New estimates have recently been given by Lipunov et al. (1995) 
and K. Postnov at this conference. 

The sequence of events accompanying the merging of two NSs was con- 
sidered by Eichler et al. (1989), Haensel et al. (1991), and Narayan et al. 
(1992). Besides being a model for GRBs (Paczyhski's earlier suggestion), 
Eichler et al. speculated that merging NSs would be sources of r-process 
elements, neutrino (z/) bursts and gravitational waves. The generation of 
the fireball is a consequence of prolific v production in binary coalescence 
through the reaction vv — >e + e~ <-> 27, as had been highlighted earlier for 
collapse events by Goodman et al. (1987). The v and v originate from pp 
and pn bremsstrahlung during the contorting tidal heating events just pre- 
ceding binary coalescence - precisely the same reactions which so quickly 
cool a NS just after birth. 

The efficiency of this reaction to produce pairs is estimated at 

(Piran et al., 1992), where Q v is the total energy in v emitted with average 
energy E v in a region with size scale Rq during injection timescale Sti, f is 
a geometrical factor, and a Dv ~ 10 _44 (-B/10 MeV) 2 cm 2 is the cross section 
for the reaction vv — >e + e~. 

Of the 10 53 -10 54 ergs of energy released in binary NS coalescence, some 
2 X 10 53 ergs appear in the form of v and ~ 8 X 10 52 ergs in the form of 
gravitational waves (Harding, 1994; see also Narayan et al., 1992). With 
£ ~ 10~ 3 , this implies that ~ 10 50 ergs is injected as fireball energy on the 
very short (~ ms) event during which the final coalescence occurs and the 
bulk of the energy transfer is accomplished, in large measure of course by 
gravitational wave damping. Numerical simulations (M. Ruffert, this con- 
ference), however, are finding smaller conversion efficiencies. If these low 
efficiencies are correct, then beaming of the pair fireball may be required to 
produce the observed 7 bursts, with the attendant required increase in p. 

7. Not-so-simple Fireballs 

During the coalescence event, tidal heating not only generates a huge lumi- 
nosity in v, but also a photon-driven mass loss when the radiation pressure 
ejects a wind. The tidal heating rate 

R* K D 1 qt 6 5 K J 

(Meszaros and Rees, 1992a), where D = 10 6 -D6 cm is the separation radius, 
iZ* is the stellar radius, to is the orbital time scale, and q is a correction fac- 
tor. An ejected mass AM ~ 1O _3 M is calculated. Only > fewxlO _9 M 
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is, however, sufficient to impede 7-ray production, because then the bulk of 
the fireball energy goes into the kinetic energy of the baryons (Shemi and 
Piran, 1990) and the burst becomes radiatively inefficient. 

Baryon contamination can be avoided by anisotropic baryon expulsion. 
Meszaros and Rees (1992b) suggested that a baryon-free zone would exist 
above and below the contact point of two coalescing NSs, and in this region 
the vv production would eject a pair jet. They also suggested that gravi- 
tational bending of v and v trajectories would produce a baryon-free zone 
in the side opposite the NS during a NS-BH merger. Mochkovitch et al. 
(1993) argue that a NS would be completely deformed into a torus dur- 
ing its merger with a BH, so that a natural axis of symmetry is formed 
along which vv annihilation occurs and a relativistically expanding e + e~ 
wind is expelled. More detailed calculations by Mochkovitch et al. (1995) 
for prescribed geometries of the ^-production region showed that the pair 
luminosity per unit solid angle L„ p ^ e + e - ~ 10 50 ergs sr _1 . The associated 
mass loss for different geometries could produce a situation where 77 decreas- 
es or the baryon loading increases with increasing v temperature in the v 
production region. 

More detailed studies of fireball evolution have been treated both ana- 
lytically (Piran, 1994; Katz, 1994a) and numerically (Meszaros et al., 1993). 
The important point is that the fireball goes through both a radiation- and 
a matter-dominated phase, as in the early universe. In its late evolution it 
also goes through a SN-like deceleration phase when the fireball sweeps up 
a sufficient quantity of matter. The fireball deceleration radius occurs at 

(Rees and Meszaros, 1992; Meszaros and Rees, 1993a), with a deceleration 
time scale 

t d ~ O.IQ^V 1 /^- 8/3 s. (10) 

Here T t = 10 3 r3 is the terminal Lorentz factor of the fireball before entering 
the deceleration phase, and n is the density of the surrounding interstellar 
medium. Some 80% of the mass is concentrated in a thin shell whose lab 
frame width remains constant until the shell becomes transparent. A reverse 
shock into the fireball accompanies the blast wave as it moves into the ISM. 
The reconversion of the kinetic energy of the baryons into radiation after 
interaction with the external medium helps relieve the problem of baryon 
contamination. 

Spectral modelling results (Meszaros and Rees, 1993b) for blast wave 
models of GRBs so far rely on synchrotron and synchrotron self-Compton 
processes, and the results are therefore very sensitive to the assumed mag- 
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netic field B. Meszaros et al. (1994) consider spectral production in three 
scenarios. For the frozen-in field model, the ejecta is assumed to have an 
entrained i?-field which is some fraction of equipartition with the total ener- 
gy density of the fireball. Dynamo amplification of B is produced by a reverse 
shock. For the turbulent model, turbulent amplification of the B is produced 
both in the blast wave as well as by the reverse shock in the ejecta, but only 
in the blast wave for the piston model. Model results show extremely weak 
emission at all energies below the 7-ray regime, but no explanation is given 
concerning the rough uniformity of the vF v peaks in GRBs near 1 MeV. 
Brainerd (1994) attributes this uniformity to Compton scattering by a sur- 
rounding optically thick medium with column density Nh ~ 10 25 cm -2 , 
which might occur if the GRB sources are buried in molecular clouds in the 
cores of galaxies. 

The fine scale structure in GRB spectra could be due to Ray leigh- Taylor 
convective instabilities in the expanding blast wave (Waxman and Piran, 
1994), overtaking shocks from unsteady outflow during the event trigger- 
ing a GRB (Rees and Meszaros, 1994), or even to patchy structure of the 
surrounding ISM. The delayed GeV emission observed in the 17 Feb 1994 
GRB (Hurley et al., 1994) has led to a rethinking concerning the contri- 
bution of different phases of fireball evolution to the observed GRB time 
history (Meszaros and Rees 1994). Perhaps the main MeV portion of the 
burst is the accelerating phase, whereas the long-duration high energy tail 
occurs when the blast wave is decelerated by the surrounding medium. Katz 
(1994b) suggests that the delayed high-energy radiation is due to the inter- 
action of the blast wave with a dense cloud and the production of 7 rays 
through secondary production reactions, notably pp — > 7r° — > 2j. 

8. Other Models 

Models where GRBs are analogs of AGNs (McBreen et al., 1993; Dermer 
and Schlickeiser, 1994; Shaviv and Dar, 1995) are favored in view of the 
similarities between the two cases, such as the peaking of the vF v spectrum 
in the 7-ray regime, the rapidly varying 7-ray emission, and the strong flaring 
behavior of blazar AGNs. Roland et al. (1994) have developed this model 
most fully and calculate complex time profiles due to the perturbation of 
the jet by a Solar mass NS or BH accreting onto a 10 3 - 1O 6 M BH. The 
differences and problems may, however, be greater than the similarities. 

The vF v spectral peaks for AGNs ranges over at least three decades in 
energy compared to only ~ 1 decade for GRBs. The compactnesses differ by 
~ 8 orders of magnitude, even for the most luminous and rapidly variable 
blazars such as PKS 0528+134. A catastrophic event, such as tidal disrup- 
tion of a star (Carter, 1992), could produce the flaring GRB behavior from 
a massive accreting black hole. But except for some carefully contrived geo- 
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metrical constructions, the minimum variability time scale corresponds to 
the dynamical time scale associated with the Schwarzschild radius of the 
BH. And if a BH with mass £ 1O 3 M is needed to agree with millisecond 
variability observed in GRBs, then this model reduces to the NS-BH merger 
scenario. Also, no AGNs have been found in GRB error boxes. 

Usov (1992; 1994) has proposed a strong field millisecond pulsar model 
for GRBs. During the collapse of a white dwarf to a NS, it is argued that 
B — > 10 15 G from flux freezing and the period P — > 1 ms from angular 
momentum conservation. For such a system, the magnetic dipole luminosity 
can exceed 10 51 ergs s _1 and the gravitational quadrupole luminosity can 
exceed 10 55 ergs s _1 . Because the energy release is so rapid, a pair fireball 
is formed, leading to a situation essentially equivalent to the NS merger 
scenarios. Collapse events leading to a pair or a Poynting-dominated MHD 
wind have also been considered by Woosley (1993) and Thompson (1994). 

9. Neutrinos from Burst Hell 

In GRB models with pion-decay (e.g., Paczyhski & Xu, 1994), the v energies 
E v range from an MeV to at least a few GeV, with flavor content in the 
ratio VeWpWr = 1:2:0. In models with emission from superconducting cosmic 
strings (e.g. Plaga, 1994), i/ e : 1/^:1/,- = 1:1:1, with E v up to ~ 10 TeV. Obser- 
vation of z/-flavor ratios could discriminate between these models. A mea- 
surement of the v to 7 flux ratio (= x) a ^ so provides important information 
about the bursting sources. Depending on optical depths, the v luminosity 
can greatly exceed that in 7 rays (x 3> 1 is possible for hadronic models 
of AGNs). On the other hand, if the 7 rays have a purely electromagnetic 
origin, x would be near zero. Weak experimental limits from underground 
experiments already exist for the v flux associated with GRBs (Miller et al., 
1994; Fukuda et al., 1994; Becker-Szendy et al., 1995), which translate into 
a limit on x ~ a few Xl0 3 (the precise limit is spectrum-dependent). 

Next year, the Super Kamiokande detector will improve the sensitivity 
for events below a few GeV by ~ 10. The effective volumes of the Baikal, 
AMANDA, NESTOR, and DUMAND ice/water instruments now under con- 
struction will be ~ 25 to 100 times that of the underground instruments, 
if the GRB v spectrum extends to E^ max > tens of GeV. The estimated 
v counting rate for each detector is ~ 10 _4 x per burst, assuming an E~ 2 
spectral shape, 1 photon (> 1 MeV) per cm 2 per GRB (Fishman et al. 1994), 
and a photon spectrum reaching E max ~ 1 TeV. With a GRB detection rate 
of ~ 1 per day, the expected number of correlated, detected z/^s per year 
is ~ 10 _2 x- Input assumptions are imprecise, so this rate could either be 
easily detectable or beyond experimental reach. 

If the energy threshold for the neutrino telescope is <C 1 TeV, then the 
expected number of v counts is ~ 5 X 10 _3 x per GRB or ~ 2x per year. 
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A possible high rate of multiple events from a single GRB is also possi- 
ble: the nearest 0.1% of GRBs might produce the spectacular signature of 
50 muons once every two years if x ~ 100. With these rates, the nature 
of the source and properties of neutrinos (masses, lifetimes, charge, speed) 
could be determined, and the weak equivalence principle could be tested by 
comparing Shapiro delays of photons, neutrinos, and anti-neutrinos passing 
the nucleus of our Galaxy. The large distance scale of GRBs allows stud- 
ies of long oscillation lengths, with possible determination of tiny v masses. 
Because the oscillation phase from z £ 0.5 sources is sensitive to the cosmo- 
logical model, v detection from GRBs could also test standard cosmology 
and further constrain Ho and Qo (see Weiler et al., 1995 for details). 

10. Finally 

The coalescing NS model remains the favorite cosmological model for GRBs 
because the sources are known to exist, the right amount of energy is 
involved, and the expected merger rates are in agreement with the required 
GRB rates. Absence of recurrence is evidently in accord with the data, 
although lensing events and triple star systems could provide escape hatch- 
es should GRB recurrence be demonstrated. Neutrino telescopes can probe 
the heart of the explosion and discriminate between models. But only the 
identification of counterparts at other wavelength ranges will conclusively 
solve the burst puzzle. 
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